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The control of biofilm formation is a challenging goal that has not been reached yet in many 
aspects. One is the role of van der Waals forces and another the importance of mutual interac- 
tions between the adsorbing and the adsorbed biomolecules ('critical crowding'). Here, a combined 
exeperimental and theoretical approach is presented that fundamentally probes both aspects. On 
three model proteins, lysozyme, Q-amylase and bovine serum albumin (BSA), the adsorption kinet- 
ics is studied. Composite substrates are used enabling a separation of the short- and the long-range 
forces. Though usually neglected, experimental evidence is given for the influence of van der Waals 
forces on the protein adsorption as revealed by in situ ellipsometry. The three proteins were chosen 
for their different conformational stability in order to investigate the influence of conformational 
changes on the adsorption kinetics. Monte Carlo simulations are used to develop a model for these 
experimental results by assuming an internal degree of freedom to represent conformational changes. 
The simulations also provide data on the distribution of adsorption sites. By in situ atomic force 
microscopy we can also test this distribution experimentally which opens the possibility to e.g. 
investigate the interactions between adsorbed proteins. 
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Proteins are found to be involved in inter- 
actions with solid surfaces in numerous nat- 
ural events. Due to protein adsorption and 
subsequent bacterial adherence, biofilms are 
formed at interfaces between solid substrates 
and liquids containing biomacromolecules. 
Biofilm formation as well as protein adsorp- 
tion are molecular assembly processes occur- 
ring at interfaces [l]. Moreover, these pro- 
cesses are collective phenomena, since surface 
coverage, adsorption kinetics, conformational 
arrangement etc. are depending on the na- 
ture of the molecular neighborhood. The lack 
of effective control over this process has be- 
come a bottleneck impeding the development 
of many biotechnologies, examples includ- 
ing biosensors, enzyme immobilization in bio- 
catalysis, antibody attachment in immunoas- 
says, biomaterial development, and tissue en- 
gineering. During the past decades, substan- 
tial progress has been made in understand- 
ing some mechanisms of protein adsorption 
PHS] . Adsorption is the net result of various 
types of interactions which depend on the na- 
ture of the protein and of the substrate, as 
well as on the surrounding aqueous solution. 
Several driving forces have been identified 
for the protein adsorption process, includ- 
ing dehydration of the protein and substrate 
surface, redistribution of charged groups in 
the interfacial film, and structural rearrange- 
ments in the biomolecules. Regarding the role 
of substrate characteristics, both surface to- 
pography and surface chemistry have been 
shown to affect protein adsorption, bacteria 



and cell adhesion [H[7j. However, since most 
topographical variations are accompanied by 
chemical heterogeneities, separating both ef- 
fects is difficult and a differentiation between 
the impacts of short- and long-range forces 
is one step towards gaining control of biofilm 
formation and has only recently been probed 
[HI [S] ■ It is the aim of this paper to compile 
the results in order to stress that the sub- 
surface composition is a parameter of utmost 
importance. 

The separation of these short- and long- 
range interactions is possible using compos- 
ite samples: By variation of the silicon ox- 
ide thickness on top of Si wafers, the van der 
Waals interactions can be tuned pUHlS) . as 
described in sectionll] The adsorption kinetics 
of BSA and a-amylase are sensitive to these 
variations of the long-range potential. Note 
that the van der Waals interaction depends 
on the geometry of the interacting objects 
[14] : For two atoms, the van der Waals po- 
tential decreases like 1/r^, where r denotes 
the distance. However, for a spherical parti- 
cle interacting with a planar surface the de- 
crease is much slower i.e. like 1/r. Hence, the 
strength of the van der Waals forces cannot 
be neglected. 

Studies concentrating on the kinetic be- 
havior during adsorption, desorption, and 
exchange processes revealed that there is a 
certain hierarchy or intermolccular arrange- 
ment of proteins in the adsorption layers (e.g. 
Schmidt et al. [15]). Additionally, reconfig- 
uration or reorientation of the adsorbed pro- 



tein can take place [111 HZ]- However, fur- 
ther information is needed to describe the 
structure of the adsorption layer in detail. 
A manifold of proteins have actually been 
detected to undergo conformational changes 
or reorientations when adsorbed on a vari- 
ety of hydrophilic, hydrophobic, charged, and 
uncharged surfaces (e.g. Keselowsky et al. 
E], Daly et al. [ig, Ishiguro et al. [Tfj . 
Baujard-Lamotte et al. [E]). These con- 
formational changes and reorientations were 
interpreted to be a result of mutual interac- 
tions between the adsorbed protein molecules 
("critical crowding"), as well as interactions 
between the proteins and between the pro- 
teins and the surfaces. Hence, monitoring the 
conformational change in an adsorbed protein 
as a function of the adsorption amount and 
time should be useful for obtaining informa- 
tion about the intermolecular interactions in 
the adsorption layer. 



The direct detection of conformational 
changes is a challenging task. Most meth- 
ods are restricted to indirect observations. 
When monitoring the adsorption kinetics, dif- 
ferent kinetic behaviour is found for the 'rigid' 
lysozyme [20] than for a-amylase and BSA 
[8l|9]. The latter is considered as a 'soft' pro- 
tein pD^E5] . This difference in the adsorption 
kinetics is therefore ascribed to dissimilar ten- 
dencies to undergo conformational changes 
due to adsorption. Monte Carlo simulations 
including an internal degree of freedom to 
incorporate conformational changes qualita- 
tively reproduce the experimental findings. 



Most methods, like ellipsometry or scat- 
tering techniques provide information that 
is obtained by averaging over some /xm^ . 
Therefore, these methods cannot provide ex- 
plicit information about single proteins in- 
volved in the adsorption process. A tech- 
nique that overcomes these shortcomings is 
atomic force microscopy (AFM). By in situ 
AFM in aqueous environment, the spatial dis- 
tribution of individual adsorbed proteins can 



be determined IV and can therefore give new 
insights into the evolution of protein films. 
The surface mobility of the single proteins 
can be tested as well as their tendency to 
form clusters or to adsorb rather independent 
from each other, as described in the scenario 
of random sequential adsorption (RS A) which 
is presented in section HI A 



I. MATERIALS 

The proteins a-amylase (product no. 10092), lysozyme 
(product no. 62971) and bovine serum albumin (prod- 
uct no. A3059) were purchased from Sigma Aldrich, 
Steinheim, Germany. BSA (66kDa)L24j and a-amylase 
(58kDa)[25] have about the same size but different iso- 
electric points (pi) at pH 4.7 [IB] and 6.5 [5S], respec- 
tively. For bovine serum albumin (BSA) the ability of 
showing different conformations is well known [2UU23J . 
Lysozyme has a mass ^27j of 14kDa and its isoelectric 
point [55] is at pH 11. 

As substrates for ellipsometry measurements, sili- 
con surfaces with different silicon dioxide thicknesses 
were used: natural (2nm, Wacker Siltronic AG, 
Burghausen, Germany) and thermally grown (192 nm, 
Silchem, Freiberg, Germany) dioxide layers, denoted as 
thin and thick oxide layers in the following. Before 
usage, the wafers must be cleaned to remove residues 
from the polishing procedure (mostly hydrocarbons). 
Therefore, the wafers were immersed into fresh 1:1 
H2SO4(conc.)/H2O2(30%) solution for 30 min. Then the 
acids were removed by rinsing the wafers in hot deion- 
ized water, three times 20 min each. All samples were 
prepared in a clean room environment (class 100), result- 
ing in a water contact angle of 0° directly after cleaning. 
Nevertheless, upon usage, the silicon wafers showed ad- 
vancing water contact angles of less than 10° (see Ta- 
ble IT] for details) and are therefore termed 'hydrophilic'. 
To gain a second set of wafers with a 'hydrophobic' 
surface (8 > 90°), a self-assembled monolayer (SAM) 
of silane molecules with a CH3 tailgroup (octadecyl- 
trichlorosilane, OTS, Sigma Aldrich, Germany) was ap- 
plied to the wafers following standard procedures [2^150] . 
The water contact angles on hydrophobized wafers were 
lager than 110° (Table IT]) with a hysteresis smaller than 
5°. Contact angle hysteresis describes the difference be- 
tween advancing Qa and receding 0^ contact angles. By 
additionally measuring the contact angles for glycerol 
and 1-bromonathalene, the surface energies as well as 
the Lifshitz-van der Waals and the Lewis acid-base com- 
ponents can be determined [31] as listed in Table lllj 

Zeta-(<^)-potential measurements [32] on wafers with 
thin and thick oxide layers with and without the OTS 
monolayer were performed for pH values ranging from 
2.0 to 7.5 (see Fig. fl]). Obviously, the hydrophilic sur- 
faces carry a higher negative charge than the hydrophobic 
wafers. The oxide layer thickness, however, does not play 
any role and zeta potential values are identical within the 
experimental error. 

The thicknesses of the silicon oxide and OTS layer were 
determined by ellipsometry, revealing an oxide thickness 
of 2 nm for the thin and 192 nm for the thick oxide layers. 
The thickness of the OTS-SAM is found to be 2.6 nm. 
The surface roughness of the four wafer types was deter- 
mined by atomic force microscopy (AFM). The measured 
root mean square (rms) roughness of an (1 /xm)^ scan 
area was below 0.2 nm for all wafer types (see Table IT] for 
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TABLE I. Silicon wafer characterization results. 
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TABLE II. Surface energy 7 and their Lifshitz-van der Waals 7^^ and Lewis acid-base 7 ^ components of the substrates, as 
determined by contact angle measurements, the results of which are listed in Table [I] 
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FIG. 1. ^-potential of hydrophobic and hydrophilic silicon 
wafers with thin and thick oxide thickness as function of pH. 
For the ^-potential, differences in oxide layer thickness are ir- 
relevant. Note that adsorption experiments in this study have 
been performed at pH 7.0 and 4.75 (vertical dashed lines). 



details). 

Keeping all other parameters (pH, ionic strength, tem- 
perature and protein concentration) constant, the four 
wafer types allow for the separation of effects due to the 
long-range and the short-range part of the surface poten- 
tial [TIM5]: Hydrophobization by the OTS-SAM alters 
the short-range part as can be seen in the contact angle 
and the C-potential. For a substantial alteration of the 
long-ranged van der Waals potential, the OTS-SAM is 
too thin. This holds for all molecular-sized layers. The 
variation of the oxide layer thickness, however, does not 
have an influence on the short-range potential, since the 
surface chemistry is identical. Hence contact angle as 
well as ^-potential measurements show no difference on 
thin or thick oxide layers. Vice versa, choosing differ- 
ent subsurface wafer compositions (in our case thin and 
thick oxide layers) enables an independent variation of 



the long-range part, i. e. the van der Waals interaction 
to probe its influence on the formation of protein films. 
The above described four types of Si wafers were used 
for probing protein film formation by ellipsometry. For 
the in situ AFM experiments shown in the following, 
mica was chosen as a substrate. Mica is a sheet sili- 
cate that can easily be cleaved. When freshly cleaved, 
mica exhibits an atomically smooth hydrophilic surface. 
Therefore, mica is frequently used for protein adsorption 
studies via AFM P3n38j. 



II. ADSORPTION KINETICS 

A. Adsorption studies by ellipsometry 

Ellipsometry [33 HO] was used to determine the ad- 
sorbed amount of proteins by in situ monitoring the 
adsorption process. The single wavelength ellipsometer 
(EP'^, Nanofilm, Gottingen, Germany) was operated in 
PCSA (Polarisator-Compensator-Sample-Analyzer) con- 
figuration at a wavelength of 532 nm. The ellipsometric 
angles ^ and A were recorded via the nulling ellipsome- 
try principle with a sampling rate of 1.5 to 6 per minute. 
This rate was sufficient to monitor the formation of pro- 
tein layers in situ. 

To determine the physical properties of the refiecting 
surface from ^ and A, a model has to be applied |39) . 
For single wavelength measurements, one has to assume 
the layers to have a constant height to be able to deter- 
mine the refractive index and thickness of all layers. For 
layers with a thickness below 5nm, it is not possible to 
distinguish between a change in refractive index and in 
film thickness [H]. Therefore, de Feijter's method is ap- 
plied to determine the adsorbed amount F as a function 
of the refractive index Uf of the protein film with a fixed 
film thickness df 



df 



r>'f 



dn/dc 



(1) 



where Ua denotes the refractive index of the ambient and 
dn/dc the increment of the refractive index of the solution 
due to the increase of molecule concentration. The refrac- 
tive index n{c) is assumed to be a linear function[321 US] 
with a fixed gradient of 0.183 cm'^/g. 

The measurements were carried out in a temperature 
controlled closed fluid cell made of Teflon. A connection 
to a flow system enables a continuous exchange of fluid 
with constant flow rates and the injection of protein so- 
lution via a sample injector (Rheodyne Manual Sample 
Injector). 

The proteins were soluted in a 10 mM phosphate buffer 
solution (pH 7.0, ionic strength 20 mM). The fluid sys- 
tem was filled with the same buffer and the buffer was 
run through the fluid cell at a constant flow rate. Af- 
ter reaching a thermal equilibrium (either at 37° C or at 
room temperature) associated with a constant baseline, 
the adsorption measurement was started by injecting the 
protein. 



B. Results and discussion 

The adsorption kinetics are shown in Fig. [2] They 
display a typical feature of adsorption kinetics: after 
a fast initial adsorption, the adsorption rate decreases 
until a plateau value of adsorbed amount is reached. 
For lysozyme (Fig. 2 (a)), the decrease of adsorption 
rate is continuous. This Langmuir-like adsorption ki- 
netics is well described in the literature [5]. The ab- 
solute adsorbed amount depends mostly on the chem- 
istry of the surface. The adsorbed amount of lysozyme 
is higher on hydrophilic samples than on hydrophobic 
ones. This can be explained by considering the posi- 
tive net charge^ of lysozyme (-1- 8 e) at pH 7 as well as 
the negative charge of the silicon wafers. The hydrophilic 
wafers (-103 mV) carry a higher negative charge than the 
hydrophobic ones (-75 mV), see Fig. fll Thus, a stronger 
attractive Coulomb interaction promotes the adsorption 
on hydrophilic wafers. For lysozyme, we observe no influ- 
ence of the thickness of the oxide layer on the adsorption 
kinetics. 

For a-amylase (Fig. 2(b)), the adsorption on thick ox- 
ide wafers also exhibits a Langmuir-like behavior [HI [9] . 
On thin silicon oxide, however, the transition between 
fast initial adsorption and saturation is not continuous. 
A regime of constant adsorption rate separates the regime 
of initial adsorption from the saturation. The final ad- 
sorbed amount depends only weakly on the oxide thick- 
ness but on the hydrophobicity of the surface. In the 
case of a-amylase, the adsorbed amount is higher on hy- 
drophobic wafers than on hydrophilic ones due to the hy- 
drophobic effect [33]. The higher attraction towards the 
surface originates from minimizing the contact of water 
with hydrophobic side groups of the protein and with the 
hydrophobic surface. As the pH of the solution is close 
to the isoelectric point of a-amylase (pi 6.9), the protein 
carries no net charge. Therefore, the Coulomb interac- 



tion just plays a minor role for adsorption affinity. 

The differences in adsorption behavior between 
lysozme and a-amylase are believed to arise from dif- 
ferent conformational stabilities [Ij. Changes of confor- 
mation due to adsorption are frequently described in the 
literature [5D|. As compared to the small and compact 
lysozyme, a-amylase can be ascribed as a 'soft' protein 
[2]. Therefore, major conformational changes are more 
likely for a-amylase than for lysozyme. 

The third protein in this study, BSA, is known for 
its ability to change conformation [201 - 125] e.g. due to 
changes of the pH of the surrounding solution. As shown 
in Fig. [31 BSA exhibits a similar adsorption kinetics as a- 
amylase |i3j . On wafers with thin silicon oxide layer (here 
only experiments on hydrophobic wafers are shown) a 
regime of constant adsorption rate is observed. When 
decreasing the temperature from physiological 37° C to 
room temperature, the regime of constant adsorption 
rate is elongated. The initial adsorption, however, is only 
slightly slowed down. Although BSA carries a negative 
net charge at pH 7 (in contrast to a-amylase) , it displays 
the same two types of adsorption kinetics as a-amylase. 
Repeating the experiments with acetate buffer at the pi 
of BSA reveals the same results. Therefore, this effect 
cannot be explained by charge effects and Coulomb inter- 
actions. The appearance of a linear regime must rather 
depend on the variation of the long-range part of the in- 
teraction potential i.e. the van der Waals interaction as 
well as the conformational stability of the protein. 

At this point it has to be stressed that numerous ex- 
perimental studies focussing on thin film formation or 
dynamics usually do not give details on the subsurface 
composition of the samples in the study. The subsur- 
face composition therefore is one to date mostly uncon- 
trolled sample property. It is very likely that many di- 
verse results about biofilm formation of different groups 
world-wide may be traced back to using different compos- 
ite samples (still exhibiting the identical surface chem- 
istry). We therefore pledge for a detailed characteriza- 
tion of samples, including e.g. a possible layering. For 
the widely used Si wafer, for instance, it is inevitable to 
characterize its oxide laver[45|. 

Since an experimental study of possible conformational 
changes of individual proteins during the adsorption pro- 
cess is too challenging to date, colloidal Monte Carlo 
(MC) simulations were launched using an effective parti- 
cle model. Here, the sensitivity of the adsorption kinetics 
to conformational (or even simpler: geometrical) changes 
of adsorbing objects shall be probed and the results shall 
then be compared to the experimental findings in order 
to be able to propose a suitable model. At first sight, 
a colloidal approach seems too rough to give reliable re- 
sults for protein adsorption. Yet after checking possible 
alternatives - which will be done in the following - it shall 
become clear that - to match the time scale of the exper- 
iment, to be able to include conformational changes and 
to use to date computer power - colloidal MC simulations 
are the method of choice. 
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FIG. 2. Adsorption kinetics of (a) lysozyme and (b) amylase on the four different types of surfaces. 
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FIG. 3. BSA adsorption kinetics on hydrophobized silicon wafers with (a) thick and (b) thin oxide layers. Measurements at 
room temperature (RT) and 37.5° C are shown. A regime of constant adsorption rate can only be observed on wafers with thin 
oxide layers and the intermediate linear regime is elongated at room temperature. 



III. MONTE CARLO APPROACH 

The numerical investigation of biofilm adsorption is a 
challenging task since the number of degrees of freedom is 
huge if the detailed structure of the proteins is taken into 
account. In addition to this, proteins undergo structural 
changes at the surfaces which implies that the relaxation 
time for the adsorption of a single protein at the sur- 
face is out of scope for present simulation approaches. 
Therefore, a drastic reduction of the degrees of freedom 
is mandatory in order to simulate the adsorption of a pro- 
tein film. In this section, we will briefly discuss the scope 
of different numerical approaches to protein dynamics. 

The most fundamental methods are based on ex- 
plicit quantum mechanical calculations. Incorporating 
the distribution of electrons they are able to predict 



the forming of atomic bonds, which is indispensable for 
studying chemical reactions of proteins (see ref . |46) for 
an overview). Ab initio methods approximately solve 
Schrodinger's equation. Being able to provide very ac- 
curate results the calculations are very time consuming. 
To date, only a two-digit number of atoms can be sim- 
ulated using this ab initio approach [47l l48]. Therefore, 
procedures have been developed [49H53] simplifying the 
unessential parts of the system using a semiempirical ap- 
proach while describing the important parts using quan- 
tum mechanical methods. In these cases, the solvent is 
partly implicit. 

Much larger systems can be considered by using frag- 
ment methods. Retaining the accuracy [54], these divide 
the whole system into several subunits which can be pro- 
cessed individually enabling the use of parallel comput- 



ers. Differing in the underlying method, several fragment 
methods are known, e.g. FMO (fragment molecular or- 
bital) method [551 lill) SFM (systematic fragmentation 
method) TT-W^ and EFP (effective fragment potential) 
method [60j .61, • Fragmentation methods have the advan- 
tage that they are nearly fully quantum mechanical in 
nature, with classical approximations often used for long- 
range interactions. They can be used for the simulation 
of small proteins: Ishida analyzed [^ the role of a polar 
protein environment in a special enzymatic process using 
FMO. Ishikawa et al determined 63] the molecular inter- 
actions between the prion protein and GN8 (a potential 
curative agent for prion diseases) based on the fragment 
molecular orbital method, too. Due to the limitation of 
the system size and time scale, however, they can't be 
used in order to study protein adsorption involving con- 
formational changes. 

To overcome those limitations, quantum effects are of- 
ten neglected. This is also the case for so-called full 
atomistic molecular dynamic simulations, where atoms 
are the smallest building blocks. By using these classical 
approaches, the simulation of peptides and small pro- 
teins becomes feasible. This includes forecasts of struc- 
tural changes of single proteins and forecasts of possible 
transition paths between two states [Ml |65] . It also in- 
cludes the simulation of changes of conformation, such as 
simpler ones that involve only loop motions. Those take 
place on a nanosecond time scale, while large changes 
that involve the rearrangement of domains may occur on 
the microsecond to millisecond time scale [66,. These 
larger changes are still not accessible using an atomistic 
level of detail. 

The time evolution of a system can be obtained via 
Molecular Dynamics, Brownian Dynamics or local Monte 
Carlo simulations, but the latter two have a limited abil- 
ity of simulating the solvent, as they can only do it im- 
plicitly. Often this is done in MD, too, as means of 
speeding up a simulation, retaining the atomistic pre- 
cision of the protein model. Whether this is reason- 
able for incidents like protein folding or adsorption is an 
open quest ion [57]. Some results suggest that it is good 
enough [5S], at least for folding. Another comparison of 
the influence of the level of detail of the solvent on the 
adsorption of peptides on hydrophobic surfaces is done 
by Sun and coworkers [551 [70] who find, that implicit sol- 
vent models yield very different results as compared to 
the explicit solvent model which was used as a reference. 
Also the different implicit solvent models predict differ- 
ent binding energies with varying level of agreement I69j . 

Further simplification for speeding up the simulations 
is done by combining multiple atoms to a single pseudo 
atom. This coarsed system needs to be described using 
modified force fields or interaction energies, respectively. 
For a fixed system size, the accessible time scales are 
two to three orders of magnitude higher than for a full 
all-atom description [71]. so in principle they allow the 
simulation of large proteins. So far, however, the cur- 
rently available generic force fields don't support struc- 



tural changes in the tertiary structure |71[. Therefore, 
reliable simulations can be obtained only if the structure 
of the protein is maintained. An overview of those so 
called united atom or. coarse-graining methods is given 
inref. [71 [73]. 

Classical MD methods are based on the integration of 
Newton's equations. The numerical effort is very high, 
typical time scales are from 10 ns [74] up to some 10 /iS 
[571 [55] , depending on the system size, the hardware be- 
ing used and accuracy (ex- or implicit water) of the calcu- 
lations. Raffaini et al for example were able to carry out 
simulations of the adsorption of single protein domains 
with a- [75] and /3-sheets [75]. Later on they analyzed [77] 
the successive adsorption of protein fragments. The ad- 
sorption of a complete but small protein, lysozyme, in- 
cluding its conformational changes was simulated by Ku- 
biak and Mulheran [7S]. Details of the large scale MD 
simulations can be found in ref. [75] . 

Omitting the time evolution completely, greatly sim- 
plifying the solvent and the intermolecular interactions 
to the point of completely stiff proteins, interaction ener- 
gies used to be calculated since two decades [50]. There 
are more recent publications using this approach[5Tl [52] 
proposing docking paths and preferred orientations, but 
conformational changes cannot be predicted by this pro- 
cedure. 

The range of applications for full atomistic as well as 
coarse-grained MD simulations indicate that these ap- 
proaches cannot be used in order to simulate the adsorp- 
tion of a full biofilm. 

Speeding up the simulations by describing the interac- 
tion with the solvent molecules stochastically using the 
Langevin-equation results in a Brownian Dynamics sim- 
ulation. Random forces with zero mean value mimic col- 
lisions of solvent molecules with the beads of a protein. 
Ravichandran et al[83] were able to show the possibility 
of adsorption of proteins with a positive net charge to a 
positively charged surfaces due to inhomogeneous charge 
distribution. Since then, however, the number of pub- 
lications using Brownian Dynamics at atomistic level of 
detail is declined in favor of Molecular Dynamics. 

Monte Carlo simulations provide another stochastic 
ansatz. They can be used both for time evolution of a sys- 
tem (when using local moves solely and detailed balance 
to reach a stationary state) or transition paths between 
two states. Simulations of coarse-grained lattice proteins 
were published some time ago: Anderson et al[55 stud- 
ied the adsorption of a model protein to an oil/water in- 
terface and Castells and coworkers '5F analyzed surface 
induced conformational changes. Being more widespread 
and potentially more accurate, today's researchers tend 
to prefer the usage of MD simulations as a workhorse 
when generating the time evolution at atomistic level. 
Nevertheless, due to their flexibility, Monte Carlo simula- 
tions offer an alternative method: Knowing an initial and 
a final state of a protein, it is efficiently possible to sample 
the transition path between those two states [5^ [551 187] . 
Sadly, the position of atoms in a protein/surface complex 



can't be extracted experimentally [31 IHH]. Not knowing 
the adsorbed state beforehand, this method cannot be 
used for protein adsorption. 

A more detailed review of the methods used for the 
simulation of conformal changes can be found in ref. [SS] ■ 



A. Collodial approaches 



Coulomb interactions in a double layer, van der Waals in- 
teractions and additionally steric repulsion. Those publi- 
cations analyze adsorption phenomena in general, in the 
latter, the parameters where chosen to mimic the adsorp- 
tion of small globular proteins such as lysozymc. Spon- 
taneous ordering depending on the ionic strength of the 
solvent is observed, too |1011ll04j . 



Coarsening the system under consideration even more, 
colloidal length scales of a few nm are reached. This 
results in a substantial increase of both spacial and tem- 
poral time scale, the latter up to the order of seconds 
to hours, depending on the level of complexity. In this 
context, the solvent is only implicitly taken into account 
and proteins typically are modelled as particles of high 
symmetry, e.g. spheres, allowing many thousands of pro- 
teins to be considered. The downside is the neglect of 
all their structural details including the distribution of 
charges and other details of the interactions which are 
relevant for very short distances. Yet, this is a necessary 
compromise in order to reach time scales that are experi- 
mentally relevant for the adsorption of a protein-biofilni. 

A simple but very efficient approach to protein ad- 
sorption is the so-called Random Sequential Adsorp- 
tion (RSA), which typically is a 2D model describing 
the parking problem omitting the particle's trajectory 
to the surface. In the simplest case, at every time 
step a particle tries to adsorb irreversibly to a random 
site by chance. The probability is determined using 
the Metropolis-Hastings algorithm in conjunction with 
a hard core potential. So the surface coverage only de- 
pends on the volume exclusion. The jamming coverage of 
this problem is described by Hinrichsen et al^U\ and the 
maximum theoretical coverage by Toth |91l . There are 
several enhancements made to this simple model. The 
adsorption and desorption of particles in combination 
with a non-spherical geometry and a change of confor- 
mation was considered by van Tassel et al 92]. Also more 
complex potentials [93H95] . e.g. for soft-RSA[96 l l97 ] are 
taken into consideration. Furthermore, the effects of 
quasi-3D models using convex particles were studied by 
van Tassel and coworkers }9 2 [ l98] . Lavalle et al[99] an- 
alyzed an RSA variant, in which adsorbed particles are 
immobilized not before a certain amount of time, allow- 
ing them to find a more suitable spot minimizing the free 
energy of the adsorbed particles. Finally, there is a more 
recent attempt to model the reversible adsorption of bi- 
nary mixtures [TOT resulting in a simple formula for the 
prediction of steady state coverage of mixtures depending 
on the steady state coverage of the pure components. In 
any case, however, the comparison to experiments is lim- 
ited, as fewer experimental parameters are realizable and 
the particle's trajectory to the surface is not accounted 
for. 

On the colloidal level, often Brownian Dynamics sim- 
ulations are used. In this case, interactions are typ- 
ically within the DLVO framework |10H - fTU5] including 



B. Simulation of a biofilm using colloidal models 
with internal degrees of freedom 




FIG. 4. Being adsorbed to a surface, particles may change 
their conformation by chance from a compact to an ex- 
panded state and vice versa using probabilities given by the 
Metropolis-Hastings algorithm. 

Using local Monte Carlo simulations, it is easily possi- 
ble to include internal degrees of freedom into the model. 
We use such an internal degree of freedom to model a re- 
versible change of conformation of a protein from a com- 
pact to an expanded state on adsorption to the surface 
as shown in Fig. l4J The expanded state optimizes the 
surface interaction at the expense of the particle/particle 
interaction whereas the compact state minimizes the cov- 
ered surface. Using this particle model we observe [HI H] 
a considerable deviation from simple, Langmuir-like ad- 
sorption kinetics (cf. Fig. [5|. 

In this specific model, DLVO interactions where uti- 
lized and adjusted using experimentally known parame- 
ters to mimic amylase adsorbing to silicon wafers. How- 
ever, due to the simplicity of the protein model many 
results are at least on qualitative level valid for other 
types of proteins as well. 

The size of the simulation box is tiny compared to the 
volume of the experimental fluid cell. Applying a canon- 
ical ensemble for the whole simulation box, the drift of 
the protein concentration in the solvent would be largely 
overestimated. Therefore it was necessary to use two sub- 
volumes (cf. Fig. pi). The upper part was to mimic an 
infinite reservoir of particles using a grand canonical en- 
semble, whereas the lower part being in contact with the 
surface had to be simulated using a constant number of 
particles, which was only changed by particles entering 
or leaving the box by means of diffusion. 

Particles adsorbing to the surface are directed to their 
adsorption site by the interaction with other particles 
already adsorbed. This is a major difference compared 
to RSA, where possible adsorption sites are chosen by 
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FIG. 5. The adsorption kinetics of simulations with (black) and without (red) change of conformation. For lower surface 
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chance. Fig. [7] shows two snapshots of the energy land- 
scape a particle moves in when approaching the surface. 
In the second plot, there is enough free space, so that a 
particle can adsorb to the surface. As a consequence, the 
energy barrier evoked by the particles already adsorbed 
is lower than in the first plot. Passing this barrier, a par- 
ticle is trapped in an attractive potential, as illustrated 
in Fig. |8] Therefore it has to following the iso line of the 
minimum and thus, the particle is directed to the adsorp- 
tion site. Actually FigJTJb) evolved from [7 a): Adsorbing 
to the surface, a particle moved it's neighbours a bit, so 



that it could fit to the adsorption site. The particle was 
removed again for the plot to ensure a better comparison 
of the energy landscape. Typically the comparison with 
the experiment is achived via the measured kinetics. It is 
suggested by experimental results that the deviation from 
Langmuir kinetics is due to long-ranged interactions with 
the silicon substrate below the oxide layer. Additionally, 
simulations indicate a similar deviation caused by con- 
formational changes which are induced by long-ranged 
van der Waals interactions. So, combining experimental 
and theoretical results, we suppose that those van der 
Waals-induced conformational changes are the cause for 
the experimentally observed kinetics. Yet, as the same 
theoretical curves could possibly be generated in different 
ways, comparing the distribution of the adsorbed parti- 
cles in both simulation and experiment would be desir- 
able. Using ellipsometry, their positions are not avail- 
able. Therefore, statistical properties like the pair corre- 
lation or Euler characteristics cannot be extracted from 
these experimental measurements. To overcome this de- 
ficiency in situ AFM measurements can be performed, as 
presented in the next section. 



IV. IN SITU AFM IMAGING 

Atomic force microscopy (AFM) was used to image the 
surface topography of the substrate as well as the ad- 
sorbed proteins in situ. For the first experiments, mica 
was used as substrate due to its simple handling and its 
smooth surface. The measurements were performed in 

TM 

TappingMode to minimize the influence of the scan- 
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FIG. 7. The energy landscape of a particle approaching the 
surface for two different points in time. The isolines mark 0.2 
steps, the white area has infinite potential energy. 




FIG. 8. The energy landscape of a particle approaching the 
surface. This is a plot of the potential energy for different 
lateral distances as shown in Fig. iTlb). 



ning tip on the proteins. A triangular cantilever (OTR 
8 or SNL, Veeco, Santa Barbara) was oscillated at a fre- 
quency of about 9 kHz. The scanning rate was between 
0.5 to 1.5 scan lines per second depending on the scan- 
ning size (ranging from (0.5 /im)^ to (3/im)^). The lat- 
eral resolution of the scans was either 256 or 512 pixels 
per line. The measurements were carried out in a closed 
fluid cell (Model MTFML, TappingModeTMFluid CeU, 
Veeco, Santa Barbara) using a MultiMode (Nanoscope 
III, Veeco, Santa Barbara). The cell was connected to 
a flow system (same set-up as in ellipsometry measure- 
ments) to allow for the exchange of buffer and protein 
solution. Experiments were performed at room temper- 
ature. Mica being frequently used for protein adsorption 
studies via AFM [55H55] was chosen as a substrate. 

The mica was mounted in the fluid cell directly after 
cleavage and consequently exposed to acetate buffer of 
pH 4.75 at a concentration of 10 mM and ionic strength 
of 8 mM. A control scan of the mica surface was taken 
before exposing the surface to protein solution. At in- 
tervals of 5 to 10 minutes, the flow was stopped and to- 
pographic scans of the surface were taken (see Fig. 10 1. 



To ensure the absence of scan artifacts, as often reported 
|105H108| . the scan area was enlarged several times to 
check for possible scan damage due to the tip. Under 
the experimental conditions mentioned above (pH, ionic 
strength) , no influence of scanning on the distribution of 
the adsorbed proteins was observed. However, increas- 
ing the ionic strength of the solution lead to signiflcant 
disturbance of the protein pattern upon scanning. 

By in situ scanning, we avoid artifacts arising from the 
use of linker molecules like glutaraldehyde which typi- 
cally lead to network-like structures (see Fig. l9lc)) or 
drying. Both procedures are frequently used in most 
AFM studies on protein adsorption to enable scanning 
in air [351 [371 [Ml L109J . 

The single objects imaged on the mica can be assigned 
to single proteins by their dimensions. The bright objects 
in Fig. [To] have an ellipsoidal shape with a long axis of 
22±7 nm and a short axis of 15±5 nm (without tip decon- 
volution: nominal tip radius 15 nm). These dimensions 
are in good agreement with the dimensions reported in 
the literature |2] (14 x 4 x 4nm'^) indicating that single 
BSA molecules are displayed on the surface. The height 
of the proteins ranges between 0.7 and 1.6 nm depending 

TM 

on the scanning parameters in TappingMode AFM. 
Therefore it is not a reliable measure for such small and 
soft objects. 

For the measurements on mica, no surface mobility was 
detected. The interaction between the negativly charged 
mica and the BSA is believed to hinder the mobility of 
the proteins. This attractive forces could be screened by 
solutions of higher ionic strenght. However, reducing the 
attraction between the proteins and the surface leads to 
fragile scanning conditions, as mentioned above. 

The proteins were not found to preferentially form 
clusters on the surface but rather form a random dis- 
tribution [llOj. To objectively characterize the distri- 
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FIG. 9. The effect of finker molecules for a stabilization of a protein film: (a) Initial configuration in the MC simulations 
using the standard form of interactions, (b) Stationary configuration in the computer simulation after adding a short ranged 
contribution to the particle-particle interactions to mimic the effect of a linker, (c) Experimental view of an amylase film 
stabilized by glutaraldehyde. For a simple comparison with the simulations, a threshold has been applied to the hight scale 
such that protein is black and substrate is white. 




FIG. 10. Series of consecutive AFM images {lfimx2fim) 
of BSA (20 /il of 0.1 mM protein solution in acetate buffer 
of pH 4.75 and ionic strength 1=8 mM at room temperature) 
adsorbed onto mica taken at the same area (slight drift) under 
stopped flow conditions. The bright objects represent single 
proteins. 



bution of these proteins on the surface, a two dimen- 
sional Minkowski analysis will be used. The Minkowski 
measures (in two dimensions: covered area, boundary 
length and Euler characteristics) are quantities that in- 
tegrate over all n-point correlation functions. This leads 
to an excellent signal to noise ratio and allows for the 
structural analysis also for a small number of samples. 
The Minkowski analysis has successfully been applied to 
characterize point patterns in biology jlllj . astronomy 
[112) and dewetting processes [113) . The knowledge of 
the statistics of adsorption sites will open an opportu- 



nity to determine the interaction potential between the 
adsorbed proteins and their surface mobility by compar- 
ison of experimental results and simulations. 



V. SUMMARY AND CONCLUSIONS 

The experimental and theoretical characterization of 
protein biofilms is a challenging task. In this work we 
used a number of different methods in order to char- 
acterize the adsorption kinetics of proteins at different 
kinds of surfaces. The experimental results showed that 
both the long- and the short-ranged interactions influ- 
ence strongly the adsorption kinetics of the proteins. By 
means of ellipsometry, it has been shown that flexible 
proteins as e.g. amylase and BSA show an intermediate 
linear regime upon adsorption on Si wafers with thin ox- 
ide layer. Complementary computer simulations, which 
describe the proteins as single spherical particles with 
an internal degree of freedom indicate that this kind of 
kinetics is the result of conformational changes induced 
by density-density fluctuations. This scenario has been 
tested experimentally by using different kinds of proteins 
and substrates. Yet the adsorption of more rigid pro- 
teins - like lysozyme - led to a Langmuir-like adsorption 
kinetics on all Si wafer types. 

A future challenge is to in situ study the adsorption 
site statistics as the surface coverage increases, which en- 
ables a direct comparison of experimental and theoretical 
results. A proof-of-concept is given by in situ AFM scans 
of BSA on mica. Future experiments will explore regimes 
of different surface mobility of the adsorbing molecules. 

To sum up, three points are important for future 
biofilm studies: i) Van der Waals forces must be taken 
into account (therefore a detailed characterization of 
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samples is mandatory), ii) conformational changes upon 
adsorption can be triggered by particle-particle and by 
particle-surface interactions (where long-ranged forces 
have to be taken into account, too) and iii) colloidal sim- 



ulations are currently the only method of choice for de- 
scribing the kinetics of these systems, which spans over 
minutes to hours. 
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